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I.  INTRODUCTION 

The  first  part  of  this  report  will  present  flowcharts, 
sample  outputs,  listings  and  descriptions  of  the  computer 
programs  used  to  produce  2-  and  4-  hour  climatic  condition- 
al probabilities  in  the  format  shown  in  Figs.  12  and  13  on 
pages  85  and  86  . Limitations  imposed  by  Saint  Louis 
University's  CDC-3300  computer  system  required  us  to  de- 
vise five  separate  programs  for  producing  the  climatic 
tables  of  this  research.  In  so  doing  each  successive  pro- 
gram uses  as  its  input  the  magnetic  tape  output  from  a 
previous  one.  Each  program  also  produces  an  audit  listing 
to  provide  an  intermediate  check  on  the  products  of  the 
system. 

Program  names  are  usually  acronyms  conceived  by  com- 
bining words  which  indicate  a general  purpose  for  the  pro- 
gram. Table  1 indicates  the  program  acronyms  used  in  this 
research  and  the  meaning  of  each. 

PGM  NAME  MEANING 

1 EXTRACTS  - Extract  data 

2 COMPUNCD  - Compute  unconditional  probabilities 

3 SMTHUNCD  - Smooth  unconditional  probabilities 

4 COMPCOND  - Compute  conditional  probabilities 

5 PRINTALL  - Print  all  conditional  probabilities 

Table  1.  Program  Acronyms 

In  this  documentation  an  attempt  has  been  made  to 
standardize  nomenclatures.  For  example,  the  array,  XPROB, 
designates  either  the  taped  output  of  program,  COMPUNCD, 
or  the  taped  input  to  the  next  successive  program,  SMTHUNCD. 
To  designate  whether  the  output  from  the  two-or  four-hour 
version  of  a prior  program  is  being  read  by  a subsequent 
one,  a suffix  of  either  2HR  or  4HR  is  appended  to  the  vari- 
able name  to  define  the  time  interval  involved.  Ceiling 
and  visibility  data  are  distinguished  in  each  case  by  the 
suffix  CIG  or  VIS. 

The  following  information  concerning  the  CDC-3300 
unique  routines  BUFFER  IN,  UNITSTF , EOFCKF  and  DECODE  are 
provided.  One  or  more  of  these  routines  may  be  found  at 
various  places  in  each  program. 

1)  BUFFER  IN  (i,p)  (a,b) 

i Logical  tape  unit  being  read, 
p Direction  and  mode  of  read. 

-1- 


J 


ialu:.  ■ 
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p = 0 Forward  read,  BCD  mode 

1 Forward  read.  Binary  mode 

2 Reverse  read,  BCD  mode 

3 Reverse  read.  Binary  mode 

a First  variable  of  the  block  to  be  transmitted 

b Last  variable  of  the  block  to  be  transmitted 

The  BUFFER  IN  statement  transmits  one  physical 
record  of  information  in  mode  p from  file  i to 
storage  locations  a through  b. 

2)  UNITSTF  (i) 

i  Indicate-;  the  Logical  Unit. 

For  the  function  UNITSTF,  the  value  returned 
is  as  follows: 

1 Buffer  operation  not  complete 

2 Buffer  operation  complete  and  no 
errors  occurred. 

3 Buffer  operation  complete,  but  an 
end-of-file  has  been  sensed. 

4 Buffer  operation  complete,  but  a 
parity  error  has  occurred. 

3)  EOFCKF  (i) 

i Indicates  the  Logical  Unit. 

EOFCKF  checks  the  status  of  the  previous  I/O 
request  on  logical  unit  i to  determine  if  an 
end-of-file  was  encountered.  The  value  re- 
turned is  as  follows: 

1 An  end-of-file  was  encountered  on 
the  last  read  operation 

2 No  end-of-file  was  encountered. 

NOTE:  The  Computed  GO  TO  statement  provides 

a convenient  method  for  checking  the 
value  returned  by  UNITSTF  and  EOFCKF. 

4)  DECODE  (c,n,v)  list 

The  DECODE  statement  converts  and  edits  informa- 
tion from  records  consisting  of  c consecutive  BCD 
characters  (starting  at  address  v)  according  to 
format  list  n and  stores  it  in  the  I/O  list 
indicated. 


Each  program  used  in  the  procedure  is  discussed  on 
the  following  pages  where  a purpose  and  description, 
flowchart,  program  listing  and  sample  output  listing  are 
given.  Finally,  a flowchart  of  the  entire  system  is 
shown  and  discussed. 

Portions  of  the  system  which  result  from  the  Saint 
Louis  University's  CDC-3300  computer  limitations  are 
cited  as  well  as  uniqueness  of  notation.  For  example, 
the  special  character  (=£)  as  listed  in  the  program 
FORMAT  statements  is  the  CDC-3300  printer  character  for 
the  character  ( ' ) . 

PROGRAM  EXTRACTS 

The  first  requirement  is  to  select  those  elements 
needed  to  produce  the  Climatic  Tables  from  the  hourly 
history  tapes.  This  program  is  designed  to  access  the 
data  base  and  select  the  following  items. 


ELEMENT 


CONTENT 


1 Initial  Observation  Time 

2 Initial  Wind  Category 

3 Initial  Dew-Point  Spread  Category 

4 Initial  Ceiling  Category 

5 Final  Ceiling  Category 

6 Initial  Visibility  Category 

7 Final  Visibility  Category 

Table  2.  Individual  observation  elements 

required  for  Stochastic  process 


The  word  'initial'  denotes  the  observation  at  the 
time  of  the  forecast.  The  word  'final'  designates  the 
observation  two-or  four-hours  later  depending  upon  the 
length  of  forecast.  The  data  card  variable  IHOUR  is  used 
to  establish  the  final  time  (See  Table  17,  page  88  ). 

EXTRACTS  is  designed  to  access  the  ETAC  TDF-14  data 
base.  With  minor  modifications  other  data  bases  such  as 
the  ETAC  DATA  SAVE  or  ARPA  DATA  BASE  may  be  used.  Only 
hourly  observations  should  be  used  since  no  final  cate- 
gories could  be  obtained  from  special  observations. 

Each  of  the  seven  fields  contain  coded  values  which 
are  used  to  correspond  to  actual  category  values  as 
follows : 


1)  Initial  Time:  The  coded  values,  0 to  23,  repre- 

sent the  actual  hour  of  the  observation  expressed 
as  Local  Standard  Time.  Care  must  be  taken  when 
using  this  value  as  an  index  to  insure  the  values 
used  are  1 to  24. 


2)  Initial  Wind  Category:  The  values,  1 to  9 , are 

used  to  express  the  wind  fields  as  follows: 


CODE  DIRECTION 


1 

2 

3 

4 

5 

6 

7 

8 
9 


0-3  KTS 
327-11 
12-56 
57-101 
102-146 
147-191 
192-236 
237-281 
282-326 


Table  3.  Wind  directions  and 
corresponding  wind 
category  codes. 


3)  Initial  Dew-Point  Spread:  The  17  dew-point  spread 

categories  are  separated  as  follows: 


CODE 

D.P.S. 

CODE 

D.P.S. 

1 

0 

10 

11-12 

2 

1 

11 

13-14 

3 

2 

12 

15-16 

4 

3 

13 

17-18 

5 

4 

14 

19-21 

6 

5 

15 

22-24 

7 

6 

16 

25-30 

8 

7-8 

17 

+ 30 

Table  4.  Dew-point  spreads  and 

corresponding  dew-point 
category  codes. 


The  card  input  variable  ITEMP  is  used  to  indi- 
cate whether  the  values  of  Table  4 are  Centigrade 
or  Fahrenheit. 
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4)  Initial  and  Final  Ceiling  Categories:  A total  of  30 

coded  values  are  used  to  express  the  ceiling  cate- 
gories as  follows: 


CODE 

CIG  HT 

CODE 

CIG  HT 

1 

0 

ft 

16 

2000 

ft 

2 

100 

ft 

17 

2200 

ft 

3 

200 

ft 

18 

2400 

ft 

4 

300 

ft 

19 

2600 

ft 

5 

400 

ft 

20 

2800 

ft 

6 

500 

ft 

21 

3000 

ft 

7 

600 

ft 

22 

3500 

ft 

8 

700 

ft 

23 

4000 

ft 

9 

800 

ft 

24 

5000 

ft 

10 

900 

ft 

25 

6000 

ft 

11 

1000 

ft 

26 

8000 

ft 

12 

1200 

ft 

27 

10000 

ft 

13 

1400 

ft 

28 

14000 

ft 

14 

1600 

ft 

29 

20000 

ft 

15 

1800 

ft 

30 

+30000 

ft 

Table  5.  Ceiling  height  values  and  corres- 
ponding ceiling  category  codes. 

5)  Initial  and  Final  Visibility  Categories:  The  30 

values  used  to  represent  the  Initial  and  Final 
Visibility  categories  are  given  in  Table  6. 


CODE 

VISBY 

CODE 

VISBY 

1 

0 mi 

16 

1 

5/  8 mi 

2 

1/16  mi 

17 

1 

3/4  mi 

3 

1/8  mi 

18 

2 mi 

4 

3/16  mi 

19 

2 

1/4  mi 

5 

1/4  mi 

20 

2 

1/2  mi 

6 

5/16  mi 

21 

3 mi 

7 

3/8  mi 

22 

3 

1/2  mi 

8 

1/2  mi 

23 

4 mi 

9 

5/  8 mi 

24 

5 mi 

10 

3/4  mi 

25 

6 mi 

11 

1 mi 

26 

7 mi 

12 

1 1/8  mi 

27 

9 mi 

13 

1 1/4  mi 

28 

14  mi 

14 

1 3/8  mi 

29 

25  mi 

15 

1 1/2  mi 

30 

+ 30  mi 

Table  6.  Ceiling  height  values  and  corres- 
ponding visibility  category  codes. 


Places  for  1000  observations  containing  each  of  the  7 
elements  are  allowed.  When  an  output  array  is  full,  all  ob- 
servations are  written  to  tape.  Since  the  last  array  may 
not  contain  1000  observations  a counter,  NUM,  is  output 
to  indicate  how  many  observations  are  contained  in  each  array 
output.  (The  limiting  factor  of  1000  is  used  so  as  not  to 
generate  an  array  too  large  for  the  CDC-3300.  For  larger 
computers  a larger  limiting  factor  should  be  used). 

The  program  requires  the  use  of  two  separate  subroutines. 
A discussion  of  each  follows. 

1)  LTRNR  (Letter-Number):  This  subroutine  is  used 

to  decode  the  over-punched  fields  of  temperature 
and  dew-point.  All  TDF-14  temperature  fields 
contain  three  digits.  Each  field  is  stored  such 
that  the  units  digit  has  a plus  or  minus  sign 
punched  over  the  digit  to  indicate  whether  the 
temperature  is  positive  or  negative.  Thus  the 
fields  must  be  broken  down  into  a two  digit 
numeric  field  and  a one  digit  alpha  field.  This 
subroutine  separates  the  alpha  digit  from  the 
overpunch  and  computes  the  desired  correct  tempera- 
ture value. 

2)  RDTAPE  (Read  Tape):  This  subroutine  is  used  to 

access  the  data  tapes.  As  previously  stated 
EXTRACTS  is  designed  to  use  the  ETAC  TDF-14  data 
base.  If  other  data  bases  are  used,  this  subroutine 
would  require  modification  to  handle  the  new  data 
base  format.  Two  points  of  caution  are  to  be  noted. 
First,  the  main  program  is  designed  such  that  one 
call  to  RDTAPE  returns  one  full  day's  data.  Any 
revision  must  take  this  into  account.  Second,  the 
TDF-14  data  base  is  such  that  all  hours  are 
accounted  for.  (Missing  observations  are  zero 
filled.)  Should  another  data  base,  such  as  the  ETAC 
DATA  SAVE  be  used,  a provision  to  handle  missing 
observations  must  be  included. 

The  following  tape  unit  assignments  are  used  by  this 
program. 

UNIT  CONTENTS 

1 Output 

2 Input  Data  Base  Tape  1 

3 Input  Data  Base  Tape  2 

4 Input  Data  Base  Tape  3 


Table  7.  Program  EXTRACTS  Input/Output 
Tape  Unit  assignments. 


A 


-7- 

The  card  input  variable  IEOF  is  used  to  indicate  the  total 
number  of  input  data  tapes  to  be  used  (See  Table  17, 
page  88  ) . 

The  following  11  pages  contain  the  flowchart,  program 
listing  and  audit  listing  for  this  program.  The  audit 
listing  indicates  the  DTG  of  the  first  and  last  observa- 
tions on  each  input  tape  and  total  count  of  observations 
output. 


Pig.  1.  "’lowchart 
for  Program  EXTRACT; 


i 


see  PROGRAM  DOCUMENTATION  FOR  DFSCRJPTION  OF  PROGRAM  FLO* 
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III.  PROGRAM  COMPUNCD 

Once  the  required  elements  have  been  extracted  from 
the  observations  and  the  appropriate  category  values 
assigned  (See  Tables  3,  4,  5 and  6)  the  frequency  of 
occurrences  for  the  various  categories  can  be  computed. 

The  observations  are  stratified  by  wind  category, 
hour,  dew-point  spread  and  ceiling/ visibility  category 
in  this  order.  Because  of  computer  limitations  a search 
is  made  of  the  entire  input  tape  for  those  observations 
which  belong  to  a single  wind  category.  A 3-dimensional 
array  (hour,  dew-point  spread,  ceiling/ visibility ) con- 
taining the  frequency  of  occurrence  values  is  computed. 
After  one  category  is  computed  the  input  tape  is  re- 
wound and  the  next  category  is  considered.  (Larger  com- 
puters which  may  use  4-dimensional  arrays  need  only  pass 
the  tape  once.) 

After  one  wind  category  is  complete  the  total 
number  of  cases  of  each  dew-point  spread  category  by  hour 
is  computed.  Values  are  cumulated  for  each  of  the  30 
ceiling/ visibility  categories  to  insure  that  all  totals 
are  increasing  with  height/distance.  The  frequencies  in 
each  ceiling/ visibility  category  indicate  cumulative 
occurrences  at  and  below  the  category  level.  Next  the 
probabilities  are  computed  by  hour  for  each  dew-point 
spread  and  ceiling/ visibility  category.  The  frequency 
of  occurrences  for  each  dew-point  spread  category  and  the 
stratified  probabilities  for  each  given  hour  are  written 
to  tape.  After  all  24  hours  have  been  written  to  tape 
for  one  wind  sector  an  end-of-file  is  written  and  the 
program  continues  to  the  next  subset.  A single  output 
tape  contains  a total  of  nine  files. 

The  following  output  tapes  are  created  from  the 
aforementioned  input  tapes. 

TAPE 

1 

2 

3 

4 

5 

6 

7 

8 

Table  8. 


TYPE 


HR  Initial 
HR  Final 
Initial 
Final 
Initial 


HR 

HR 

HR 


HR  Final 
HR  Initial 
HR  Final 


Ceiling 

Ceiling 

Visibility 

Visibility 

Ceiling 

Ceiling 

Visibility 

Visibility 


The  eight  different  type  of 
runs  required  by  this  program 
for  the  Stochastic  process. 
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The  tape  unit  assignments  for  this  program  are  as  follows: 

UNIT  CONTENTS 

1 Input  (EXTRACTS  Output) 

2 Output 

Table  9.  Input/Output  tape  unit  assign- 
ments for  program  COMPUNCD. 


The  next  seven  pages  contain  the  flowchart,  program 
listing  and  a sample  audit  listing  for  this  program.  The 
audit  listing  indicates  the  total  number  of  observations 
by  hour  for  the  indicated  wind  category  and  dew-point 
spread. 
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Fig.  4.  Sample  output  Audit 
listing  from  program  COMPUNCD, 
Indicates  by  hour/ dew- point 
spread  the  total  frequency  of 
occurrence  for  wind  category. 
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PROGRAM  SMTHUNCD 

This  program  is  designed  to  smooth  the  unconditional 
probabilities  output  by  program  COMPUNCD.  Each  day's 
data  are  read  into  an  array  similar  to  that  used  in 
COMPUNCD  except  the  hour  dimension  is  now  26.  Each  hour's 
data  are  read  into  elements  2 through  25,  leaving  elements 
1 and  26  unfilled. 

The  frequencies  are  stored  in  a 3-dimensional  array 
with  the  x-direction  being  the  24  hours,  the  y-direction 
the  17  dew-point  spread  categories  and  the  z-direction  the 
30  ceiling/visibility  categories.  All  smoothing  is  done 
in  the  x-y  plane  for  each  of  the  30  ceiling/ visibility 
categories  of  the  z-plane.  As  a natural  result  of  the 
hour  dimension  being  cyclic  in  nature,  the  probabilities 
along  the  hour  edges  are  exchanged  to  perform  better 
smoothing  (e.g.,  the  values  for  hour  24  which  are  stored 
in  element  25  are  placed  in  element  1 and  those  for  hour 
1 which  are  stored  in  element  2 are  placed  in  element  26). 

Each  data  point  is  assigned  a weighting  factor  for 
the  x-direction  (hour)  and  the  y-direction  (dew-point 
spread).  The  dew-point  weighting  factor  is  determined  as 
follows : 

DPWGHFAC  = ( 0 . 114  3*N  + 1.0)*5, 

where  N is  the  total  observations  for  the  particular 
spread.  The  hour  weighting  factor  is 

HRWGHFAC  = 0 . 6667*(DPWGHFAC  - 1.0)  + 1.0 

Each  ceiling/ visibility  category  is  processed  in 
such  a way  that  the  probability  at  each  point  uses  eight 
surrounding  probabilities  to  obtain  a smoothed  proba- 
bility value.  As  an  illustration  of  the  smoothing  scheme, 
data  point  (2,2)  would  be  smoothed  as  follows: 

P(l,l)  + H( 1 ,2 ) *P( 1 , 2 ) + P(1 ,3)  + 

D(2 ,1)*P(2 ,1)+D(2 ,2)*H(2,2)*P(2 ,2 )+D( 2 , 3 ) *P( 2 ,3)+ 
P( 3 ,1)  + H(3,2)*P( 3,2)  + P( 3,3) 

P( 2 ,2  ) = 

1.0  + H( 1 ,2 ) + 1.0  + 

D(2  ,1)  ♦ D( 2 ,2 ) *H( 2 , 2 ) + D(2,3)  + 

1.0  + H( 3 ,2 ) + 1.0 

Where  P represents  the  probability  of  the  point,  D the 
dew-point  spread  weighting  factor  and  H the  hourly 
weighting  factor.  The  formula  shows  that  each  new  proba- 
bility is  based  upon  the  summation  of  the  probabilities 
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of  each  point  times  a weighting  factor  divided  by  the  sum 
of  the  weighting  factors.  The  weighting  of  each  corner  is 
taken  to  be  one.  The  edges  are  smoothed  using  only  the  five 
surrounding  points. 

The  program  does  a preliminary  smoothing  in  those  areas 
of  the  array  where  the  initial  probability  is  zero.  The 
initial  smoothing  cycles  a maximum  of  ten  times  or  until  all 
values  are  non-zero.  All  probabilities  are  converted  to 
exponential  form  before  entering  the  main  smoothing  scheme. 

The  main  smoothing  scheme  cycles  for  a variable  number  of 
tines  as  set  by  the  input  variable  ILIM.  (Investigation 
has  shown  that  ILIM  = 8 produces  the  optimum  smoothing. ) 

Finally,  when  all  30  ceiling/ visibility  categories  have 
been  smoothed,  the  probabilities  for  each  hour  are  written 
to  tape  along  with  the  frequencies  at  each  dew-point  spread, 
ceiling/ visibility  intersection.  Note  sample  output  in 
Fig.  7.  The  frequencies  written  are  those  which  have  been 
smoothed  and  not  the  unsmoothed  input  frequencies.  This  is 
done  so  as  to  have  the  smoothed  values  which  are  to  be  used 
to  compute  the  unconditional  probabilities  in  program  COMPCOND 
and  displayed  in  the  rubric  of  the  Climatic  Tables. 

This  program  uses  one  subroutine,  PRTDATA  (Print  Data), 
to  list  either  the  input  unsmoothed  unconditionals  (Fig.  6) 
or  the  output  smoothed  unconditionals  (Fig.  7).  The  card 
input  variable  IHOUR,  ISEASN,  ITYPE , IMODE , ITEMP  and  ISTN 
select  the  correct  heading  to  be  printed.  The  card  input 
variable  IPRT  selects  how  many  hours  of  each  day's  data  are 
to  be  printed.  A maximum  of  24  x 9 x 12  = 432  separate 
matrices  could  be  printed.  This  includes  2 matrices  (un- 
smoothed and  smoothed)  for  each  24  hours  and  9 wind  categories. 
Since  one  matrix  output  requires  2 pages,  the  total  listing 
would  contain  864  pages. 

The  following  tape  unit  assignments  are  used  by  this 
program. 


UNIT  CONTENT 

1 Input  ( COMPUNCD  Output) 

2 Output 

Table  10.  Input/Output  tape  unit  assign- 
ments for  program  SMTHUNCD. 

The  next  14  pages  contain  the  flowchart,  program  listing 
and  a sample  output  of  the  smoothed  and  unsmoothed  proba- 
bilities computed  by  this  program. 
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LOOP  THROUGH  ALL  *IND  CATEGORIES 
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RESET  NEW  PROBABILITIES  ARRAY  INTO  OLD  ARRAY  FOR  FUTURE!  SMOOTHING 


NOW  SET  ALL  PBOBAR IL I T I PS  TO  AN  EXPONENTIAL  FORM 
ANY  PROBABILITY  STILL  ZERO  IS  SET  TO  .00001. 
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Fig.  7.  Sample  output  for 
program  GKTHUNCD.  Indicates 
frequency  and  probability  ^or 
each  ceiling  category  and  Dew- 
point spread  after  smoothing. 
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The  frequencies  listed  in  the  smoothed  output  (Fig.  7) 
were  recomputed  from  the  smoothed  probabilities.  Hence 
the  frequency  values  for  specific  indices  for  the  un- 
smoothed (Fig.  6)  and  smoothed  (Fig.  7)  probabilities  may 
not  be  the  same.  Areas  of  the  display  which  are  not 
printed  have  either  frequencies  of  zero  occurrence  before 
smoothing  cr  truncate  to  zero  after  smoothing. 

V.  PROGRAM  COMPCOND 

Procedures  to  this  point  have  pertained  to  uncon- 
ditional probabilities.  Observations  were  selected  from 
a given  season  (hour,  dew-point  spread,  wind  direction, 
ceiling  and  visibility  categories)  for  both  an  initial 
and  two-or  four-hour  final  time  to  define  these  proba- 
bilities. They  were  then  subjected  to  smoothing  pro- 
cedures. We  shall  now  proceed  to  deduce  conditional 
probability  estimates  using  the  initial  and  final  (2  or 
4 hours)  unconditional  probabilities  and  the  gridded 
representations  of  the  universal  graphs  for  estimating 
conditional  climatologies,  i.e.,  the  Stochastic  model 
contained  in  the  appendix  of  the  original  report. 

The  smoothed  unconditional  probabilities  provided  by 
program,  SMTHUNCD,  are  read  from  tape  for  the  initial 
and  final  hours.  Twelve  of  the  seventeen  dew-point 
spreads  and  sixteen  of  the  thirty  initial  ceiling/ visi- 
bility categories  contained  on  these  tapes  are  required 
to  produce  the  Climatic  Tables.  (See  discussion  of  sub- 
routine SELTINT  on  page  47. ) Table  11  lists  the  required 
dew-point  spread  categories  and  Table  12  lists  the  six- 
teen initial  ceiling  and  visibility  categories  processed. 
The  five  categories  selected  from  the  array  of  Table  12 
which  are  contained  in  the  format  for  field  usage  are 
shown  in  Table  13  (See  SELTFIN). 


OUTPUT 

INPUT 

CATEGORY 

D.P.S. 

CATEGORY 

1 

0 

1 

2 

1 

2 

3 

2 

3 

4 

3 

4 

5 

4 

5 

6 

5 

6 

7 

6 

7 

8 

8 

8 

9 

10 

9 

10 

15 

12 

11 

20 

14 

12 

+ 30 

17 

Table  11.  The  17  input  dew-point  spread 

category  codes  and  the  corresponding 
output  dew-point  spread  category  codes. 
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CIG/VIS 

OUTPUT 

CIG 

CIG 

INPUT 

VIS 

VIS 

INPUT 

CATEGORY 

VALUE 

CATEGORY 

VALUE 

CATEGORY 

1 

0 

ft 

1 

0 

mi 

1 

2 

100 

ft 

2 

1/16 

mi 

2 

3 

200 

ft 

3 

1/8 

mi 

3 

4 

300 

ft 

4 

1/4 

mi 

5 

5 

400 

ft 

5 

1/2 

mi 

8 

6 

500 

ft 

6 

3/4 

mi 

10 

7 

600 

ft 

7 

1 

mi 

11 

8 

800 

ft 

10 

1 1/2 

mi 

15 

9 ■ 

1000 

ft 

11 

2 

mi 

18 

10 

1500 

ft 

13 

3 

mi 

21 

11 

2000 

ft 

16 

4 

mi 

23 

12 

2500 

ft 

18 

5 

mi 

24 

13 

3000 

ft 

21 

6 

mi 

25 

14 

5000 

ft 

24 

7 

mi 

26 

15 

10000 

ft 

27 

10 

mi 

28 

16 

NO  CIG 

30 

15 

mi 

29 

Table  12. 

The  input  ceiling/visibility  category 

codes  and 

the  corresponding  16  output  ceiling/visibility 

CIG/VIS 

category  codes. 

CIG  CIG  CIG 

VIS 

VIS 

VIS 

OUTPUT 

LETTER 

VALUE 

INPUT 

LETTER 

VALUE 

INPUT 

CATEGORY 

CODE 

(FT) 

CATEGORY 

CODE 

(MI) 

CATEGORY 

1 

A 

0- 

199  2 

J 

0 - 

7 

2 

3 

200- 

499  5 

K 

5s  - 1 

10 

3 

C 

500- 

999  10 

L 

1-2 

17 

4 

D 

1000-2999  20 

M 

2-3 

20 

5 

E 

3000-9999  26 

N 

3-6 

24 

Table  13.  The  five  final  output  ceiling  and  visibility 
category  codes. 


The  conditional  probabilities  can  be  computed  once  the 
data  for  one  hour  has  been  read  and  the  required  initial  and 
final  unconditional  probability  values  are  known.  To  do  this, 
a gridded  display  stored  in  the  computer  is  accessed  in  which 
the  initial  unconditional  probability  defines  the  ordinate, 
the  final  unconditional  probability  the  abscissa  and  the 
corresponding  conditional  probability  value  it  read-off  at 
the  intersections  of  these  coordinate  values.  As  a result 
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of  the  packing  of  these  gridded  displays  (Universal  Graphs) 
care  must  be  taken  to  insure  that  the  right  values  are  used. 
(See  appendix  for  racking  scheme.)  Because  the  graphs  only 
contain  values  from  0 to  100  in  increments  of  02  an  interpo- 
lation scheme  is  used  to  obtain  intermediate  values.  A 
check  is  made  to  insure  that  none  of  the  computed  values  are 
larger  than  100. 

In  this  interpolation  scheme,  the  conditional  probability 
values  corresponding  to  the  next  highest  and  next  lowest  even 
integer  values  (IXY1,  IXY2  , IXY3,  IXY4)  with  respect  to  a 
given  coordinate  intersection  along  the  ordinate  and  abscissa 
are  first  obtained  from  the  graphs.  The  differences  (CX, 

DX,  CY,  DY)  between  the  required  values  and  the  unconditional 
values  are  computed  and  used  to  weight  each  of  the  conditional 
probability  values  to  obtain  the  required  conditional  proba- 
bility (See  Fig.  8 for  example). 

IXY3  IXY4 

CY 

DX  CX 

DY 


IXY1  IXY2 

Fig.  8.  Values  used  in  the  interpolation 
scheme  for  computing  the  con- 
ditional probability  from  the 
Universal  Graph. 


Once  the  conditional  probabilities  for  all  dew-point 
spreads  and  initial  categories  have  been  computed  for  one  hour, 
a check  is  made  to  make  sure  that  the  cumulated  value  in  each 
final  category  is  equal  to  or  larger  than  that  of  the  pre- 
ceding category.  Next  the  cumulated  value  for  the  immediately 
lower  category  is  subtracted  to  obtain  the  actual  probability 
for  that  category.  The  probability  values  for  the  topmost 
category  (F  or  0)  is  found  by  subtracting  the  cumulated  value 
of  the  5 lower  categories  from  100.  Prior  to  writing  the 
values  to  tape,  one  last  check  is  made  to  insure  that  all 
probabilities  decrease  with  increasing  dew-point  spread. 

Finally  all  conditional  probabilities  for  one  hour  are  written 
to  tape  with  the  unconditional  probabilities  for  the  12 


dew-point  spreads  and  specific  wind  direction  (See  discussion 
on  subroutine  COMPUNCD) . The  same  procedure  is  used  for  all 
hours  and  wind  subsets. 

Program  COMPCOND  requires  four  separate  subroutines. 

Each  is  discussed  below. 

1)  SELTINT  (Select  Initial):  This  subroutine  selects 

the  12  dew-point  spreads  and  16  initial  categories 
which  comprise  the  initial  values  of  the  Climatic 
Tables  from  the  17  dew-point  spread  and  30  initial 
categories  of  the  program  SMTHUNCD.  For  dew- 
points, the  array,  KCAT,  is  used  to  select  the  12 
spread  categories.  The  array  is  composed  of  17 
elements,  one  each  for  the  17  original  dew-point 
spread  categories.  Those  dew-point  spreads  which 
are  input  and  are  to  be  used  as  an  output  category 
from  this  program  contain  a consecutive  value  to 
indicate  the  output  category  value.  Those  cate- 
gories which  are  not  to  be  output  contain  a zero 
value  and  are  skipped.  The  array,  ICAT,  is  a two 
dimensional  array  containing  30  elements.  Each 
corresponds  to  the  original  30  ceiling  and  visi- 
bility categories.  Similar  logic  is  used  for  the 
ICAT  array  to  that  just  discussed  for  KCAT.  Those 
category  elements  of  program  SMTHUNCD  which  are 

to  be  used  as  an  input  to  a later  program  contain  a 
designation  value  corresponding  to  the  category  to 
be  output.  Categories  which  are  not  to  be  used 
contain  a zero  designator.  Note  Table  12  on  page  45 
which  gives  a cross  reference  between  the  input  and 
output  category  values  and  the  corresponding  ceiling/ 
visibility  values. 

2)  SELTFIN  (Select  Final):  This  subroutine  selects 

those  dew-point  spread  and  ceiling/visibility  cate- 
gories which  are  to  be  used  as  the  final  categories 
(see  Table  13  on  page  45  of  this  report).  The  pro- 
cedure by  which  the  array,  KCAT,  is  used  to  select 
the  final  12  dew-point  spread  categories  was  pre- 
viously discussed  under  section  1 above.  The  same 
logic  is  used  in  this  subroutine.  The  array,  JCAT, 
in  this  subroutine  is  similar  to  ICAT  in  the  sub- 
routine SELTINT  differing  only  in  the  number  of  final 
categories  that  are  output.  The  five  categories  of 
JCAT  are  indicated  in  Table  13  on  page  45. 
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3)  COMPUNCD  (Compute  Unconditionals):  This  subroutine 

is  used  to  compute  1)  the  unconditional  probabili- 
ties of  the  final  categories  and  the  12  dew-point 
spreads  occurring  in  the  specific  wind  direction, 

2)  final  ceiling  probability  when  dew-point  strati- 
fication is  not  considered,  and  3)  final  ceiling 
probability  when  neither  dew-point  and  wind  strati- 
fication are  not  considered  for  the  rubric.  Note 
the  values  at  the  bottom  of  Figs.  12  and  13  on  pages 

85  and  86.  In  the  rubric  information  the  unconditional 
probabilities  for  wind  direction  are  not  stratified 
by  dew-point  spread.  They  are  computed  by  summing 
all  frequency  of  occurrence  values  for  the  required 
ceiling/ visibility  category  and  dividing  by  the 
total  occurrence  for  that  hour  and  wind  direction. 

The  computation  of  the  probabilities  in  the  ALL  WINDS 
category  is  simply  the  summation  of  the  frequency  of 
occurrence  for  each  given  temperature  dew-point 
spread  category  for  each  of  the  individual  wind  dir- 
ections divided  by  the  number  of  observations  of  that 
temperature  dew-point  spread  irrespective  of  wind 
direction.  The  ALL  WINDS  values  are  computed  subse- 
quent to  the  wind  and  temperature  dew-point  spread 
stratified  ones.  The  array  IUNPRBAW  is  written  as  the 
tenth  file  on  the  output  tape  to  output  these  values. 

4)  PRTDATA  (Print  Data):  This  subroutine  lists  the  com- 

puted conditional  probabilities  prior  to  outputting 
them  to  tape.  See  the  sample  output  for  this  program 
in  Fig.  10  on  page  65.  Similar  to  previous  programs 
the  card  input  variable,  IPRT,  is  used  to  set  the  array 
NHOUR  which  determines  how  many  hours  of  data  are  to 

be  printed.  A total  of  two  pages  are  required  to 
list  one  hour's  output  of  conditional  probabilities. 

The  tape  unit  assignments  required  by  this  program  are 
indicated  in  Table  14. 

UNIT  CONTENTS 

1 2/4  Hour  Initial  Ceiling/Visibility  data 

2 2/4  Hour  Final  Ceiling/ Visibility  data 

3 Appropriate  Universal  Graphs  tape 

4 Output  data  tape 

Table  14.  Input/output  tape  unit  assignments  re- 
quired by  program  COMPCOND. 

The  next  17  pages  contain  the  program  flowchart,  sample 
program  listing  and  sample  program  output  for  the  program 
COMPCOND. 


SEE  PROGRAM  DOCUMENTATION  FOR  DESCRIPTION  OF  PROGRAM  FLO 
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Fig.  10.  Conditional  probabilities  as  computed  by  the  program 

COMPCOKD.  Value 8 shown  are  for  the  indicated  wind 
direction  and  hour.  Left  margin  indicates  the 
initial  ceiling  category  and  final  letter  categories. 
Initial  dew-point  spread  categories  are  indicated  , 
along  the  top  margin 
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VI.  PROGRAM  PRINTALL 

This  program  is  used  to  list  the  conditional  and 
unconditional  probabilities  computed  by  program  COMPCOND 
in  the  format  of  Figs.  12  and  13.  The  two-and  four-hour 
ceiling  and  visibility  values  are  simultaneously  input 
into  the  program  for  printing. 

The  program  first  reads  through  the  nine  wind  cate- 
gories of  program  COMPCOND  to  obtain  the  unconditional 
probabilities  for  the  respective  ceiling  and  visibility 
occurrences . 

The  two-  and  four-hour  data  for  ceiling  are  processed 
first.  The  four-hour  values  must  first  be  read  since  the 
array  used  to  hold  the  unconditional  probabilities  for  the 
specific  wind  direction  and  each  of  the  12  dew-point 
spread  categories  are  used  to  hold  both  the  two-  and  four- 
hour  data.  Thus  the  valid  two-hour  data  overlay  the 
previously  read  invalid  four-hour  values  in  the  arrays. 

Once  the  values  for  the  medians  have  been  computed 
and  all  probabilities  rounded  to  one  digit  the  data  are 
listed.  The  output  is  arranged  such  that  the  conditional 
probabilities  are  printed  for  each  wind  direction  for 
all  hours  in  sequence.  For  one  hour's  display  the  ceil- 
ing data  are  listed  followed  bv  the  visibility  data.  When 
nine  wind  categories  are  considered  for  each  hour  of  the 
day,  the  total  number  of  output  pages  is  432  for  ceiling 
and  visibility  inclusive. 

A total  of  five  separate  subroutines  are  required  by 
this  program.  A discussion  of  each  follows. 

1)  MEDVALUE  (Medium  Value):  This  subroutine  is 

used  to  compute  the  MED  values  for  the  con- 
ditional probabilities  for  both  the  two-and 
four-hour  final  categories.  The  MED  value 
is  defined  as  that  level  where  the  cumulative 
probability  of  lower  ceiling  is  the  same  as 
that  of  higher  ceilings,  i.e.,  the  50th  per- 
centile level.  The  probabilities  for  each 
final  category  is  summed  until  a value  greater 
than  50  is  obtained.  Once  found  the  amount  of 
the  current  category  needed  to  make  the  50 
percentile  is  divided  by  the  probability  of 
the  category.  This  value  is  then  multiplied 
by  the  layer  increment  and  added  to  the  lower 
level  of  the  layer.  The  same  routines  apply 
when  visibility  categories  rather  than  ceilings 
are  under  consideration. 
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2)  ROUNDOFF  (Round  Off):  This  subroutine  is  used  to 

round  all  conditional  probabilities  such  that  they 
can  be  output  as  a single  digit.  All  rounded  values 
are  obtained  by  first  adding  5 and  then  dividing  by 
10.  In  this  way  an  initial  value  of  54  would  round 
to  5 and  a value  of  55  would  be  rounded  to  6.  No 
attempt  is  made  to  force  the  values  to  add  up  to  10. 
Should  a value  of  95  or  greater  be  found  for  any 
given  category,  that  category  is  rounded  down  to  9. 

3)  UNCDEWPT  (Unconditional  Dew-Point):  This  subroutine 

is  used  to  compute  the  MED  value  and  round  all  un- 
conditional probabilities  for  each  of  the  12  dew- 
point spread  categories.  The  same  procedure  as  out- 
lined in  the  subroutines  MEDVALUE  and  ROUNDOFF  is  used. 

4)  UNCWINDS  (Unconditional  Winds):  This  subroutine  is 

used  to  compute  the  MED  value  and  round  the  uncon- 
ditional probabilities  for  the  two  wind  values  listed. 
The  first  being  the  ALL  WINDS  category  and  the  second 
the  specific  wind  direction.  The  same  procedure  as 
outlined  in  the  subroutines  MEDVALUE  and  ROUNDOFF  is 
used. 


5)  PRTDATA  (Print  Data):  This  subroutine  is  used  to  list 

the  conditional,  unconditional  and  MED  values  as  com- 
puted in  the  prior  subroutines  (see  Figs.  12  and  13). 
Care  must  be  taken  in  listing  the  data  since  the  MED 
values  for  ceiling  are  integer  and  the  MED  values  for 
visibility  are  floating  point. 

This  program  requires  the  following  tape  unit  assignments. 


UNIT 


CONTENTS 


1 2 HR  Ceiling  Conditional  Values 

2 2 HR  Visibility  Conditional  Values 

3 4 HR  Ceiling  Conditional  Values 

4 4 HR  Visibility  Conditional  Values 

Table  15.  Input/Output  tape  unit  assignments 
for  program  PRINTALL. 


The  next  19  pages  contain  the  program  flowchart,  program 
listing  and  two  sample  output  listings. 
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COMPUTE 


DEW  PT 
MED  a 
ROUND 


lg.  11a.  flowchart  for  program  PRINTALL  continued. 
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Fig.  12.  A page  of  the  final  format  for  ceiling  climatologies 
at  Offutt  AFB  when  the  initial  forecast  time  is  0300  and  the 
wind  is  0-3  KTS. 


Fig.  13.  A page  of  the  final  format  for  visibility  climatologies 

at  Offutt  AFB  when  the  initial  forecast  time  is  0300  and  the  wind 
0-3  KTS. 
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VII.  SUMMARY  OF  COMPLETE  SYSTEM 

The  production  of  Climatic  Tables  such  as  those 
illustrated  in  Figs.  12  and  13  for  an  entire  season  for 
both  the  wind  and  dew-point  stratification  presented 
requires  that  the  five  previously  described  programs  be 
run  a total  of  23  times.  To  simplify  this  routine,  an 
attempt  has  been  made  to  generalize  the  data  cards  which 
are  required  by  a program.  Although  some  programs  may 
not  require  all  the  information,  each  program  reads  the 
same  formatted  data  card.  Table  17  describes  each  of 
the  input  variables  and  the  numeric  values  used  to 
indicate  the  various  parameters.  Also  indicated  are 
the  program  (see  Table  16)  which  requires  each  of  the 
parameters . 

The  chart  in  Fig.  14  depicts  the  flow  for  the  en- 
tire system.  As  can  be  seen,  each  program  is  required 
to  be  run  various  number  of  times  as  indicated  in 
Table  16. 


PGM#  NAME  RUNS 

1 EXTRACTS  2 

2 COMPUNCD  8 

3 SMTHUNCD  8 

4 COMPCOND  4 

5 PRINTALL  1 


Table  16.  Program  number,  program  name  and 
number  of  runs  required  by  each 
program. 


The  main  difference  between  the  runs  of  the  same 
program  is  in  the  input  tapes.  In  the  case  of  the  pro- 
gram COMPUNCD  the  main  difference  is  in  the  data  card 
which  is  input. 

Some  of  the  variables,  such  as  ISEASN  and  ISTN, 
will  remain  constant  for  all  runs  of  the  programs. 

Others,  such  as  IHOUR  and  ITYPE , will  vary  systematically. 
A discussion  of  each  program  and  the  card  input  values 
which  vary  follows.  (NOTE:  PGM#  in  the  following  Tables 

refers  to  the  number  found  in  the  upper  right  hand 
corner  of  the  box  containing  the  program  name  in  the 
system  flowchart. ) 
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NAME 


IEOF 


ILIM 


ISTN 


CONTENTS 

Number  of  input  tape(s) 


PGM  § 


IHOUR 

2 

for 

2-HR  increment 

1 

2 

3 

4 

4 

for 

4-HR  increment 

ISEASN 

1 

for 

Spring 

1 

2 

3 

4 

2 

for 

Summer 

3 

for 

Autumn 

4 

for 

Winter 

ITYPE 

1 

for 

Ceiling 

2 

3 

4 

2 

for 

Visibility 

I MODE 

1 

for 

Initial 

2 

3 

2 

for 

Final 

ITEMP 

1 

for 

Temperatures 

in 

°F 

1 

2 

3 

4 

2 

for 

Temperatures 

in 

°C 

IPRT 

Number  of  hours  to 

be 

printed 

per  day  (Must  be  factor  of  24) 

3 

4 

5 

Number  of  cycles  desired  in 
smoothing  (ILIM=8  is  suggested) 

Name  of  station  being  processed 
(Maximum  of  32  Characters) 


1 2 3 4 5 


Table  17. 


Data  card  input  variables,  the  numeric  values 
which  each  can  have,  and  sequential  number  of 
each  program  which  uses  the  various  input 
variables . 


1)  EXTRACTS:  This  is  the  only  program  which  uses  the 

variable  IEOF.  IEOF  is  used  to  indicate  how  many 
data  base  input  tapes  are  to  be  used.  This  program 
is  run  twice — once  to  output  the  two-hour  final 
categories  and  once  to  output  the  four-hour  final 
categories.  The  data  card  input  variable,  IHOUR , 
changes  for  the  two  runs  as  indicated  in  Table  18. 


PGM  # 

1 

2 


IHOUR 

2 

4 


Table  18.  System  flowchart  program  number  and  the 

values  IHOUR  obtains  for  each  of  the  two  runs. 
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2)  COMPUNCD:  The  system  flowchart  shows  that  this  pro- 

gram uses  the  two  output  tapes  from  EXTRACTS  and 
processes  each  of  them  a total  of  four  times.  The 
card  input  variables,  IHOUR,  ITYPE  and  IMODE  are 
used  to  indicate  whether  the  data  to  be  processed 
is  two  or  four  hour,  ceiling  or  visibility  and 
initial  or  final  data.  Table  19  depicts  the  system- 


variation 

of  the  input 

parameters 

indicated 

PGM  # 

IHOUR 

ITYPE 

IMODE 

3 

2 

1 

1 

4 

2 

1 

2 

5 

2 

2 

1 

6 

2 

2 

2 

7 

4 

1 

1 

8 

4 

1 

2 

9 

4 

2 

1 

10 

4 

2 

2 

Table  19.  System  flowchart  program  number  and 
the  values  IHOUR,  ITYPE  and  IMODE 
obtain  for  each  of  the  eight  runs  . 


3)  SMTHUNCD:  Due  to  the  paucity  of  data,  each  set  of 

the  unconditional  probabilities  must  be  smoothed  to 
make  them  sufficiently  reliable  for  producing  con- 
ditional probabilities  by  the  methods  of  this  re- 
port. Therefore,  each  of  the  eight  previous  output 
tapes  of  program  COMPUNCD  are  to  be  used  as  an  input 
to  this  program.  The  card  input  variables  IHOUR, 

ITYPE  and  IMODE  are  likewise  used  to  indicate  the 
type  of  data  being  processed.  Table  20  depicts  the 
systematic  variation  of  the  input  parameters  indicated. 


PGM  § 

IHOUR 

ITYPE 

IMODE 

11 

2 

1 

1 

12 

2 

1 

Jl. 

1 

13 

2 

2 

1 

14 

2 

2 

2 

15 

4 

1 

1 

16 

4 

1 

2 

17 

4 

2 

1 

18 

4 

2 

2 

Table  20. 


System  flowchart  program  number  and 
the  values  IHOUR,  ITYPE  and  IMODE 
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NOTE : In  this  program  and  those  to  follow  the  variable 

IPRT  indicates  the  number  of  hours  for  which  the 
data  is  to  be  printed  (e.g.  a value  of  IPRT  = 2 
would  print  HOUR:  0 and  HOUR:  12,  or  a total  of 
two  of  the  24  hours,  a value  of  IPRT  = 3 would 
print  HOUR:  0,  HOUR:  9 and  HOUR:  18,  or  a total  of 
three  of  the  24  hours). 

4)  COMPOUND:  This  program  is  run  a total  of  four  times 

to  combine  like  output  (2/4  hour-Ceiling/ Visibility) 
tapes  from  the  previous  program.  As  shown  in  Table  21, 
the  card  input  variables  IHOUR  and  ITYPE  indicate 
whether  the  data  is  two-  or  four-hour,  ceiling  or  visi- 
bility. 

PGM  # 

19 

20 
21 
22 

Table  21. 


IHOUR 


ITYPE 


2 

2 

4 

4 


1 

2 

1 

2 


System  flowchart  program  number 
and  the  values  IHOUR  and  ITYPE 
obtain  for  each  of  the  four  runs, 


5)  PRINTALL:  Finally,  the  four  output  tapes  from  the 

previous  program  are  ready  for  printing.  The  input 
variables  ISEASN  and  ISTN  indicate  on  the  Climatic 
Tables  the  season  and  station  for  which  the  data  was 
processed.  The  variable  IHOUR  indicates  the  first 
hour  of  the  two  hour  group.  In  this  way  either  2-4 
HR  or  3-6  HR  values  may  be  output  without  modifica- 
tion to  the  program.  The  variable  IPRT  should  be 
set  to  24  so  as  to  output  all  hours. 

Reference  the  numbers  inside  the  tape  symbols  on  the 
system  flowchart  (Fig.  14).  The  top  number  indicates  the  out- 
put tape  unit  assignment  for  the  previous  program  and  the 
bottom  one  indicates  the  input  tape  unit  assignment  to  the 
next  program. 

Fig.  14  on  the  next  page  contains  the  complete  system 
flowchart  to  produce  the  output  for  one  season.  Other  seasons 
are  Droduced  by  repeating  the  procedure  for  a different  value 
of  ISEASN. 


V 
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r 

VIII. 


PROGRAM  EXECUTION  TIME  REQUIREMENTS 

No  attempt  is  made  to  estimate  the  time  required 
for  each  of  the  five  programs  to  execute  on  a specific 
computer  except  for  the  Saint  Louis  University  CDC-3300. 
Only  general  comments  as  to  the  CDC-3300  execution 
speed  can  be  made  since  the  time  required  is  a function  of 
the  amount  of  data  printed  by  each  program. 

The  program  EXTRACTS  when  used  to  process  the  data 
for  Offutt  AFB  required  approximately  38  minutes  each 
for  the  one  and  a half  reels  of  input  data.  The  two 
runs  to  obtain  both  the  2 and  4 hour  final  categories 
required  about  1 hour  and  20  minutes.  Systems  which 
would  use  a data  base  other  than  the  TDF-14  should  run 
much  faster.  Since  the  TDF-14  data  tapes  are  in  Binary 
Coded  Decimal  (3CD)  the  program  execution  is  greatly 
slowed  due  to  the  conversion  from  BCD  values  to  integer 
for  every  observation. 

Each  of  the  runs  of  the  program  COMPUNCD  required 
approximately  20  minutes.  Thus  for  the  eight  runs  a 
total  of  2 hours  and  40  minutes  is  required.  Since 
larger  computers  would  need  to  read  through  the  input 
tape  only  once  the  time  should  be  decreased  by  a factor 
of  nine. 

Program  SMTHUNCD  requires  about  20  minutes  per  run 
or  a total  of  2 hours  and  40  minutes  for  the  eight  runs. 
This  estimate  is  based  on  the  time  taken  to  print  only 
one  hour’s  data  for  each  of  the  nine  wind  categories. 

Because  of  the  gross  inefficiency  of  the  program 
COMPCOND  in  using  a tape  for  the  Universal  Graphs,  this 
program  requires  one  hour  to  run  for  each  of  the  four 
required  runs.  A total  of  4 hours  is  required  for  this 
program.  Again  this  estimate  is  based  on  the  time  re- 
quired to  print  only  one  hour's  data  for  each  of  the  nine 
wind  categories.  The  use  of  a disk  to  hold  the  Universal 
Graphs  would  greatly  improve  this  program's  efficiency. 

The  program  PRINTALL  is  the  one  program  where  the 
amount  of  data  printed  greatly  determines  the  execution 
time.  To  print  all  24  hours  ceiling  and  visibility 
data  for  all  nine  wind  categories  in  their  entirety  re- 
quires about  65  minutes  computer  time.  As  stated  pre- 
viously the  total  output  consists  of  432  pages  per 
season . 


J 
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Totaling  the  above  execution  times  indicates  that 
approximately  11  3/4  hours  of  computer  time  on  a 
CDC-3300  is  required  per  station  per  season.  This  value 
might  appear  to  be  excessive  if  it  weren't  for  the  fact 
that  larger  computers  and  increased  program  efficiencies 
would  decrease  this  execution  time  substantially. 


IX.  APPENDIX 

1 ) Introduction 

This  appendix  details  the  format  of  the  magnetic 
tape  which  contains  the  matrix  values  of  the 
Universal  Graphs.  It  is  assumed  that  the  reader 
has  some  knowledge  of  magnetic  tape  storage  charac- 
teristics . 

2 ) General  Comments 

Contained  on  the  tape  are  the  matrix  values  used 
to  produce  the  120  2-4  HOUR  Universal  Climatic 
graphs  and  the  120  3-6  HOUR  Universal  Climatic 
graphs.  As  used  in  this  project  the  tape  was  written 
at  556  BPI , Binary  mode  (7-track).  The  CDC-3300  is 
a 4 byte  per  word,  6 bit  per  byte  computer.  Each 
matrix  contains  data  values  from  .00  to  1.00  in  in- 
crements of  .02.  Thus  each  unpacked  matrix  (see 
Appendix  Section  D)  contains  a total  of  51  x 51  = 2601 
data  points.  The  graphs  are  naturally  subdivided  into 
six  groups  of  five  each  as  follows: 


A 

TO 

(A) 

A 

TO 

(A-B) 

A 

TO 

(A-C) 

A 

TO 

(A-D) 

A 

TO 

(A-E) 

B 

TO 

(A) 

B 

TO 

(A-B) 

B 

TO 

(A-C) 

B 

TO 

(A-D) 

B 

TO 

(A-E) 

C 

TO 

(A) 

V* 

TO 

(A-B) 

C 

TO 

(A-C) 

C 

TO 

(A-D) 

C 

TO 

(A-E) 

D 

TO 

(A) 

D 

TO 

(A-B) 

D 

TO 

(A-C) 

D 

TO 

(A-D) 

D 

TO 

(A-E) 

r 

TO 

(A) 

£ 

TO 

(A-B) 

E 

TO 

(A-C) 

E 

TO 

(A-D) 

E 

TO 

(A-E) 

F 

TO 

(A) 

F 

TO 

(A-B) 

F 

TO 

(A-C) 

F 

TO 

(A-D) 

F 

TO 

(A-E) 

Table 

22. 

Combinations  possible 
universal  graphs. 

from  one  set 

of  the 

The  same  subdivided  groups  are  valid  for  the  visi- 
bility categories  J,  K,  L,  M,  N and  0. 

3)  File  Format 


In  order  to  calculate  the  Climatic  Conditional  Proba- 
bilities as  reproduced  in  the  sample  output  for  program 
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COMPCOND  (Fig.  10),  it  is  necessary  to  have  all  five 
matrix  of  one  category  available  in  the  computer  at  one 
time.  Thus  each  set  of  five  matrix  are  placed  on  the 
tape  such  that  a single  READ  will  transfer  all  13005 
data  point  values  into  the  computer  memory.  An  end-of- 
file  was  written  after  each  series  to  facilitate  search- 
ing for  the  required  set  of  data.  Thus  the  tape  con- 
tains six  files  for  each  of  the  eight  types  or  a total 
of  48  files.  The  eight  different  types  are  contained  on 
the  tape  in  the  following  order. 


HOUR 


TYPE 


2 

4 

2 

4 

3 

6 

3 

6 


Ceiling 

Ceiling 

Visibility 

Visibility 

Ceiling 

Ceiling 

Visibility 

Visibility 


Table  23.  The  order  of  the  individual  universal 
graphs  on  the  master  tape. 


The  program  COMPCOND  requires  that  only  the  graphs  for 
one  type  be  on  the  Universal  Graphs  input  tape. 


4 ) Date  Packing 

Because  of  the  limited  memory  capacity  of  the  Saint  Louis 
University  CDC-3300,  all  13005  matrix  values  could  not  be 
transferred  into  memory  if  one  location  was  used  for  one 
value.  Thus,  the  values  are  packed  in  such  a way  that 
three  values  occupy  the  same  location.  In  this  way  the 
resulting  4335  locations  are  within  the  memory  size  limi- 
tations of  the  computer.  As  a result  of  the  packing  the 
previous  (51  x 51  x 5)  matrix  is  transformed  into  a (17 
x 51  x 5)  matrix. 


The  data  values  are  packed  in  the  following  order.  Input 
matrix  points  (1,1,1),  (2,1,1)  and  3,1,1)  are  placed  in 
the  output  matrix  point  (1,1,1).  Matrix  values  (4,1,1), 
(5,1,1)  and  (6,1,1)  go  into  (2,1,1).  This  system  is 
carried  throughout  ail  51  rows  and  columns  of  the  five 
initial  category  matrix. 
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The  following  logic  is  used  to  place  the  three  values  into 
one.  Input  matrix  point  (3,1,1)  is  multiplied  by  10,000 
and  placed  in  the  above  indicated  output  point.  Next 
point  (2,1,1)  is  multiplied  by  100  and  added  to  the  out- 
put point.  Finally,  point  (1,1,1)  is  added  to  the  output 
point.  Again  this  scheme  is  used  throughout  all  51  rows 
and  columns. 


Sample  Listings 


The  next  four  pages  contain  a sample  listing  and  sample 
output  for  the  program  which  can  be  used  to  list  the  data 
values  from  the  Universal  Graphs  Master  Tape.  Because  of 
printer  limitations  and  the  desire  to  produce  a readily 
useable  output,  increments  of  .04  are  listed  for  the 
abscissa.  Values  of  intermediate  increments  can  be  assumed 
to  be  linear  between  listed  values. 


THIS  PROGRAM  is  DESIGNED  TO  READ  THE  UNIVERSAL  grapes  tape 
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Fig.  15.  Sample  listing  of  one  of 
Ceiling  as  produced  by  the  program 


the  Universal 
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